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Abstract 

Coupled Ising models are studied in a discrete choice theory framework, where they can be understood to represent 
interdependent choice making processes for homogeneous populations under social influence. Two different coupling 
schemes are considered. The nonlocal or group interdependence model is used to study two interrelated groups making 
the same binary choice. The local or individual interdependence model represents a single group where agents make two 
binary choices which depend on each other. For both models, phase diagrams, and their implications in socioeconomic 
contexts, are described and compared in the absence of private deterministic utilities (zero opinion fields). 

1 Introduction 

Social scientists have been interested in the formal study of discrete choice theory settings for decades. These allow for 
the consideration of a wide variety of problems, ranging from the study of social pathologies [T] to demand contexts 
[21 [3j HJ [5j HI [7] or election results [3 [9] , as well as the prevalence of different relevant social traits or opinions [5j [TUJ [TTJ 
[T21 [T51 03] ■ Traditionally, the effects of social influence were neglected. Some of the work by Schelling P2J HHJ [TT] and 
Becker [TH1IH] are exceptions, for some of which direct links to models from condensed matter can be drawn [201 I2T1 122] . 
Interactive setups similar to the ones that will be considered in this paper were introduced by Follmer |23j and Granovetter 
[23]. In the 1990s, economists such as Blume, Brock and Durlauf started developing a consistent framework in which 
to systematically study these problems taking into account social interactions and noted the deep connections to some 
models from statistical physics [251 123 HZ1 12H] • These ferromagnetic setups had also been already proposed by some of 
the first physicist advocates of building a physics like body of knowledge for the social sciences (TUJ [2H 130] • 
The traditional approach typically considers a group of agents maximising the utility function: 

U(si,hi,Ei(s),ei(si)) = u(si,hi) + S(si,Ei(s)) + ei(sj) (1.1) 

where u is the private deterministic utility, S the social deterministic utility and e, the private random utility. The choice 
Si made by the individual i can take the values -1 or 1. Each is a component of vector s, with Ei(s) representing the i th 
agent's subjective belief on the choices of the rest of the agents. The idiosyncratic willingness to adopt (hereafter IWA) 
or opinion field hi characterises personal preferences of each agent, its inherent predisposition towards the choice making 
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process. Finally, e; is an individual and choice dependent random shock. Economic demand contexts can be studied in 
these setups by interpreting IWAs as a combination 6j — p, where bi is the idiosyncratic willingness to buy (hereafter 
IWB) and p the fixed price. Demand curves, i.e., the dependence on price of the fraction of adopters fi (fraction of 
the individuals with Sj = 1), can be studied for fixed values of the rest of the parameters, as discussed, for example, in 

[21 H 0131]. 

There is a direct analogy that can be established between these settings (and their total utility) and condensed matter 
models for ferromagnetism (and their total energy). The binary choice variable Si can be understood as a binary spin 
variable and so the average choice s = Sj as the system's average magnetisation, which can be easily expressed in 

terms of the fraction of adopters fj, as s = 2/i — 1. The action of the IWA hi, and therefore of private deterministic utility 
terms, will be associated in its ferromagnetic counterpart to the action of magnetic fields. Social deterministic utility can 
be related to spin interactions. Random private utility is translated into the introduction of random noise. Assuming a 
logistic distribution for the difference of the random payoff terms P(ej(— 1) — £i(l) < z) = 1+e 1 _^ z identical for all agents, 
is equivalent to studying the ferromagnetic system using the canonical ensemble, i.e., in statistical equilibrium. Note that 
in the physical model, (3 is given the statistical mechanical usual interpretation of inverse of the temperature (3 = 1/KbT, 
were Kb is Boltzmann's constant. In discrete choice settings, socioeconomic temperature, 1/(3 in the utility discrete 
choice setting, represents statistical fluctuations that can reflect an incomplete knowledge about agent's particularities, 
individuals' probability of behaving irrationally and/or a more fundamental uncertainty concerning free will. This allows 
for an interpretation of (3 in terms of social permeability, which measures the potential response of individuals to social 
influence [TUj . 

Under particular choices for private and social deterministic utilities, the system can be described using a mean field 
Ising model with constant external field [25J [57]; that was solved exactly in the 1930s (see for example [52]). This is 
only true when considering a private deterministic utility u proportional to both the individual's choice and the IWA, 
which is considered the same for all agents (constant opinion field). Social deterministic utility will have as many 
terms as other agents with whom there is interaction. All agents are assumed to interact with all others and each of 
these terms to be proportional to the quadratic difference between the agent's opinion and its perception of the other's 
choice: S = ~Ei (^2jjti '^t ( Si — s j) 2 )- The coupling Jy measures how much agent i wants to conform to the opinion 
of agent j and thus the strength of social influence/magnetic interaction and is also assumed to be constant throughout 
the population [§. Agents have rational expectations, i.e., Ei(sj) = s. Note this means all individuals will be equally 
compelled to align their opinions with the average choice. All agents have the same (accurate) perception of what the 
average behaviour of the group is. As all agents also have identical IWAs, populations are completely homogeneous. 

These are drastic simplifications concerning most socioeconomic systems of interest but are useful to understand 
basic qualitative features emerging from social interaction. Complexity can be added and hypothesis relaxed to make 
them more realistic and parallelisms to statistical physics still remain useful. It seems particularly relevant to include 
heterogeneity that enables the characterisation of the group. This can be achieved using random fields and/or spin glass 
type models, which have been extensively studied in physics [331 1331 ESI I3S1 EZl I3B] - Work on these lines in socioeconomic 
contexts can be found in [H [5J [7J HU [551 HS1 1301 [32 133 SO] • Other interesting variations on the interacting scheme, such 
as nearest neighbours or similar [II [3j E [HI [9l [I3J [23j [26l [28j HI1 [42j dynamic interacting neighbourhoods [401144] or 
complex networks [SI 1431 145] can be used. However, the mean field approach can be argued to be a good approximation 
for many problems of interest to mimic a general tendency to conform to the norm and is more tractable analytically 

3 Dimensionless choice variables Si are considered, and so utility, social coupling, IWA and temperature are measured in the same units 
(formally equivalent to the physical scenario in natural units h = Kb = 1), which can be denned to suit the specific problem under consideration. 
They can be an attempt to quantify abstract quantities (satisfaction, wcllbeing, reputation. . . ) or measure specific gains or losses (surpluses 
or deficits) of the individuals, for example, in money. 



2 



gllZllinillliniESlEZlEHllSBISllSllinillSlllS- Relaxing rational expectations and/or reciprocity assumptions [Tl[31][46] 

are also interesting paths to explore, as is the study of finite size effects [3Ql [39l 00] . The study of this kind of setups as 
dynamical systems (which can evolve out of equilibrium) |31 [HI El HU HI 1131 EH EH EH EH EH] is currently an active area 
of research sometimes referred to as opinion dynamics (see |40j and |48j for reviews) . Similar setups can also be used to 
study other related problems such as resource allocation, hierarchical structures, coalition formation, social learning, and 
public goods games [2 [6l [20j [23l [25l [4Q1 [43l [46l [47l [4^ 

When comparing discrete choice settings subject to social influence to those without, interesting qualitatively different 
characteristics appear. Even the simple Ising mean field model with constant field already presents a rich statistical 
mechanical phenomenology [35] with interesting translations to socioeconomic language [TUJ [T2] H3 [2S1 1201 EH] ■ The 
qualitative average outcome (phase) is described by its order parameter, the average choice s. For the unpolarised or 
paramagnetic phase (s = 0), in average half of the group is deciding for and half against. In the polarised or ferromagnetic 
state (s ^ 0), there is a greater prevalence of one of the options over the other. 

Unbiased populations (with respect to a particular choice) arc those with zero IWA for all individuals. The average 
outcome will be determined only by social interaction, individuals having no inherent deterministic preference for any 
of the choices. We can for example consider this to be the case for some fashions and traditions, where there is no real 
cost or benefit for the individuals besides the social payoff provided through imitation. The unpolarised state can only 
be stable in this case, with a second order phase transition for J/3 = 1, where a system will experience a rapid but 
continuous change from an unpolarised state (for low social permeability /? or high social coupling J) to a polarised one. 
In the polarised regime (1 < J/3), the stable equilibrium is degenerate, as there are two identical maxima of the utility 
s = ±7« (minima of the free energy, physically equivalent states). The previous history of the system becomes relevant 
and microeconomic specifications of the model may not uniquely determine its macrocconomic properties. 

When we turn on the opinion fields, these critical values of the parameters for which a second order phase transition 
takes place for the unbiased case, still separate regions where social utility is relevant (spontaneous magnetisation exists 
in the physical analog) and regions where it is not. Systems in a previously unpolarised state at h = will turn into 
well behaved, single stable polarised states, with the sign of the average choice determined by that of the population's 
IWA. If the system was already in a polarised state at zero IWA, while the sign and value of the average choice will 
in general be unambiguously determined by that of the field, for small enough values of the latter, the symmetry is 
not completely broken and the degeneracy of the average choice equilibria persists. While no longer equiprobable or 
symmetric, depending on the previous state of equilibrium of the system, the final average magnetisation may not be 
aligned with the IWA. Metastable states can be understood as a collectively reinforced choice that persists even when 
the rest of conditions are ripe for a change in the sign of the average choice. In economic demand contexts, it involves 
two coexisting low and high demand states for a given price. They are associated to a first order phase transition at 
h = 0, which provides a mechanism for abrupt, discontinuous changes in the sign of the average choice and the possibility 
for hysteresis (changes in the state of the system when varying the parameters that can not be undone by reversing the 
process). 

Some attempts at explaining data referring to real social processes using this type of approach have been made 
[U [2j [5] [9] [TT] [14]. This line of research remains, however, somewhat scarce, and more efforts need to be put in this 
direction. 

In this work, besides considering social influence, two different binary choice variables, each governed by Ising mean 
field dynamics, are coupled in another ferromagnetic like interaction. These and similar setups have been studied in 
physical contexts, particularly for their interest in explaining plastic phase transitions EH EH E21 Es3 EH E5] ■ It amounts 
to the introduction of an additional contribution to the utility (jl.ip that bonuses or penalises agents for the alignment 
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of both spin variables. This is the case of interest when considering two different groups making the same choice when 
both groups are interdependent in their decision making process. Each individual is subject to social pressure to conform 
to its own group's average. Besides, agents will also tend to align (k > 0) or dis-align (k < 0) their options with the 
other group's average choice. This setup is referred to as the group interdependence or the nonlocal model (agents can 
be considered to interact with all members from its own group and the other one). Coupled spin variables are also the 
case of interest to study a single group in which each individual has to make two choices which are related, which can be 
labelled individual interdependence or local model (the coupling between both spin variables is then only through each 
individual). These systems can be in three possible phases depending on the two (coupled) order parameters: unpolarised 
or paramagnetic phase (s = t = 0), polarised or ferromagnetic phases (s ^ 0, t ^ 0) and mixed phases (one zero and one 
nonzero order parameterfl 

In the next two sections, Hamiltonians for both models are described and expressions for their free energies and 
systems of equations of state analytically derived. Section 0] describes and compares the resulting phase diagrams for the 
zero opinion fields or unbiased case. Finally, some concluding remarks are made. 



2 The group interdependence or nonlocal model 

Let us consider two groups, one made up of N s agents making binary choice Sj = ±1 and the other of N t agents making 
binary choice tj = ±1 with dynamics governed by the Hamiltonian: 

where summations over (i,j) are 1 < i < j < N s for the first term, 1 < i < j < N t for the second term and 1 < i < N s , 
1 < 3 < N t for the mixed term. Summations over i arc 1 < i < N s for the fourth term and 1 < i < N t for the last term. 

The intra- couplings J s , Jt are positive and constant and quantify the strength of social pressure within each group. 
Each individual's decision is also affected by the agent's correct perception of what the other group is doing through 
the inter- coupling k. It may either reinforce (k > 0) or discourage (k < 0) alignment of choices between individuals in 
different groups but interaction is always symmetric (exactly reciprocal). Using mean field theory, which is exact for 
infinite groups (thermodynamic limit in physical terminology), is equivalent to considering each agent as influenced by its 
correct perception of the average choice of the other group. Besides intra- and inter- group social interactions, each group 
is subject to a constant opinion field. The most natural application seems to problems where the interaction between both 
groups is of social nature. It could also refer, however, to actual gains or losses for agents unrelated to social interaction. 

This setup could be used to study problems such as the incidence of social pathologies in neighbouring districts, 
public opinion on a certain issue in two neighbouring countries, technology choice in two different related industries, 
female participation in the labour market in two neighbouring regions, etc. 

When both groups are of the same size (which is the case when considering infinitely large groups) this is the system 
studied in [52]. This work extends the results presented there using numerical methods to calculate stable solutions of 
the system in thermodynamic equilibrium. 

For N s = Nt = N, using mean field theory (on all intra- or inter-coupling terms) we can rewrite (|2.ip as: 

H = ^-^s 2 + ^y-t 2 + Nkst - (J s s + kt + h s ) ^ s i ~ ( J tt + ks + h t ) ^ U ( 2 - 2 ) 

i i 

4 This terminology convention is different from that used in 1551 , where the ferromagnetic phase is referred to as mixed phase. 
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where s = -h Si and t = ^ i ij are the average choices in each of the groups. The corresponding partition function 
for the representative canonical ensemble (Z = Tre~@ H where Tr indicates sum over all possible spin configurations) can 
be expressed 

Z = e-^f fi2 J =+f * 2 Jt+Nkst) [ 2 cogh ^ (JsS + kt + hs ^ 2 cosh (/3 (J t f + ks + h t ))] N (2.3) 

where /? is the social permeability or inverse of the temperature (which in this case accounts for statistical fluctuations) . 

Finally the system's free energy density (/ = F/N with F the free energy F = /3 _1 log(Z)) will be given in the mean 
field approximation by 

/ = 1 J sS 2 + ^J t t 2 + kst-^j In (2 cosh (/3 ( J s s + kt + h s ))) - ^ In (2 cosh (/3 (J t t + ks + ht))) (2.4) 

Stable states of the system will be those minimising the free energy. After derivation and simple algebraic manipulation, 
critical points can be found to be given by the system of equations of state: 

a(a-t&nh(P(J a a + kt + h a ))) = 
a(t-t&Dh(P(Jtt + ks + ht))) = 

with a = J s Jt — k 2 . 

There are thus two different cases. For the degenerate case (a = 0), substituting J t = J s /k 2 in the original system 
both give the same equation of state: 

J t s + M-J a tanh(P(J a s + H + h 8 ))-ktaxih(p(^t + k8 + h^ =0 (2.6) 

For the non-degenerate case (a =/= 0), on which this work focuses, critical points will be given by solutions to the 
system of equations of state: 

s = tanhp (J s s + kt + h s )} 
t = tanh[/3 (J t t + ks + h t )] 

Positive definite Hessian (evaluated at the solution of the system (j2.7|) ) is a sufficient condition for a critical point to 
be a minimum. In this case, its determinant can be written as: 

dct(H) = a (1 - pJ sls - p.Jat) + P 2 a 2 lslt (2.8) 
and the free energy's second derivative (and first diagonal element of the Hessian) as: 

-L = J S - pj s 2 ls - pk 2 lt (2.9) 
as 2 - 

where ls = cosh , w ^ s+kt+hs}] and 7t = cos ^ mj ^ t+ks+ht)r Positive and (EU) are sufficient conditions for minima. 

When zero opinion fields are considered, polarised solutions of the system of equations of state (|2.7[) will always 
appear in pairs (s, t) and (— s, — t), each pair of solutions having the same stability. There is also a symmetry relating 
solutions under a change in sign in k and in any one of the average magnetisations. Furthermore, the unpolarised state or 
paramagnetic solution s = t = is always a critical point of /. Critical regions where its stability changes can be studied 
by linearising equations (|2.7|) for \s\ 1 and \t\ I (at finite non zero temperature) when h s = h t = 0, yielding: 
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s = fi{J s s + kt) + 0{s 3 ,t 3 ,s 2 t lS t 2 ) 
t = f3{J t t + ks) + 0(s 3 1 t 3 1 s 2 t,st 2 ) 



(2.10) 



Further simplification of this system leads to the expression: 



a(3 2 -{J s + J t )p + 1 = 



(2.11) 



Solving (|2.1ip gives two values of j3 where the stability of the paramagnetic phase changes from saddle point to 
maximum / minimum: 



In the strong coupling regime (k 2 > J s Jt), there is a single physical (positive) value, where the unpolarised solution 
changes from saddle point to maximum, and so of no relevance in the discussion of phase transitions. 

In the weak coupling regime (k 2 < J s Jt), there will be two physical values of /3. The larger one (given by the + option 
in (|2.12p ) is of the same type as discussed above. The smaller one (given by the — ), is a critical point at which the 
unpolarised or paramagnetic phase will change from saddle to minimum and a second order phase transition takes place. 
For temperatures bellow the critical one, the system is in a polarised or ferromagnetic phase. This means that, even in 
the absence of personal predispositions of any type, social interaction both with individuals from within the agent's own 
group and from the other one, promotes the emergence of a tendency to polarisation (which is only complete, i.e., s = ±1 
and t = ±1 at T = 0). The actual direction of choice is not determined, only the alignment or dis- alignment between 
both average choices (there are two physically equivalent states of equilibrium) . Note that when opinion fields are turned 
on, even if the symmetry is broken and the paramagnetic solution will no longer be a critical point, this will still delimit 
regions of the phase diagram where interaction in itself has an impact on the aggregate outcome. 

3 The individual interdependence or local model 

Let us now consider a single group of N identical individuals all of which are simultaneously making two interdependent 
binary choices Sj = ±1 and tj = ±1. There will be an additional term in the utility function that penalises or bonuses 
the agent for aligning or not its two choices. Both of the decisions are subject to social influence to conform to the norm 
and to a constant opinion field. A Hamiltonian describing such a system is: 



where summations over arc 1 < i < j < N and summations over i are 1 < i < N. Note that in this case the 

interaction between both choice variables cannot be of social nature. 

Discrete choice problems of interest to the social sciences that can be explored using this model include the relation 
between teenage pregnancy and school dropout in a given population, results in simultaneous referendums or in elections in 
a bipartidist system, demand for two different software products of the same brand, motherhood and female participation 
in the labour market and so on. 

This is a similar model to that studied in |55) . where random inter-coupling k is considered. Most results presented, 
however, are for symmetric probability distributions p{k{) which seem of little interest in the socioeconomic context, as 



J a + J t ± V(J S + Jt) 2 - 4a _ J s + J t ± V(J S - J t ) 2 +4fc 2 



(2.12) 



2a 2a 




(3.1) 
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they involve populations for which half of its members have a positive interaction between both choices while for the 
other half they tend to oppose. 

Hamiltonian ()3.1[) can be rewritten: 

H= Js + Jt_Js_^ s ^ -A^Tt^) + + MO (3.2) 



Following [55] , the partition function of the representative canonical ensemble can be computed exactly for an infinitely 
large system (thermodynamic limit) yielding: 



oo />oo 



Z=^(J a J t )*e-!< J ' +J *> / 



oo J — oo 



where g(s,t) is the free energy functional such that the free energy density / = lirn/v->oo 77 = dsdtg(s,t) and 

is given by the expression: 

g = ij s s 2 + ij t i 2 - iln[2e /3fc cosh (/3 {J s s + J t t + h s + h t )) + 2e~ f3k cosh {/3 (J s s - J t t + h s - ht))] (3.4) 

Critical points of the free energy / are those of the functional g. Introducing the notation a s = tanh (/3(J s s + h s )), 
at = tanh (/3( J t t + ht)) and a k = tanh (/3k), the system's equations of state can be written as: 

a s + a t ak 

s 



1 + a s a t ak 

In this case the Hessian's determinant takes the form: 



1 + a s a t ak (o c\ 

a t + a s a k ^■ > 



det(H) = J s J t [l -— + 7— 74 (3.6) 

(1 + a s a t a k ) l (1 + a s a t a k r 

and g's second derivative with respect to s is given by: 

Pg lVhs{l-aW k ) ( , ?] 

ds 2 s [l + a s a t a k y { ' 

The zero opinion field case has the same symmetries described for the nonlocal non-degenerate zero field model. The 
paramagnetic or unpolarised solution is also always a critical point of g. Linearization of Q3.5P for \s\ 1 and \t\ <^ 1 (at 
finite nonzero temperature) when h s = h t = yields: 

s = (3J s s + /3 J t tanh(/3fc)£ + 0(s 3 , t 3 , s 2 t, st 2 ) 
t = 13 J t t + f3J s t&nh(f3k)s + 0(s 3 ,t 3 ,s 2 t,st 2 ) 

which can be simplified to 

J s J t (1 - tanh 2 (^fc)) P 2 - [J t + J s )f3 + 1 = (3.9) 

Depending on the values for the intra- and inter-couplings, there can be three qualitatively different scenarios when 
analysing equation (|3.9p . with either one, two or three different solutions in (3 (all positive and thus physically relevant). In 
all cases, only the smaller (or only) one of these (the highest temperature) will be a critical point where the paramagnetic 
or unpolarised phase becomes unstable and there is a second order phase transition taking place. 
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4 Phase diagrams for the unbiased case 



The Newton— Raphson algorithm has been used to numerically solve the equations of state for the unbiased case (h s = 
h t = 0), fixed J s = 1 in some given units (which is equivalent to measuring the rest of couplings and fields in terms of 
J s ) and different values of the rest of the parameters. Python libraries developed for the computation and analysis of 
these are available at https://github.com/anafrio/. Results are summarised through figures of the phase diagram's 
cross-sections for both models for specific values of the rest of parameters. These have been chosen for the sections to 
be representative of all the possible phenomenology. More details on the solutions found (including unstable ones) and 
dependence on the different parameters can be found in |56) . 

Unbiased situations refer to what has been loosely referred to as trends or traditions, i.e., social conventions which 
will give no particular advantage or disadvantage to the isolated individual and choice. For the group interdependence 
model, this would be the case, for example, to study the use of a particular outfit or accessory in two different socioe- 
conomic groups, or the adoption of specific vocabulary (argot, technical term) by two such groups. Unbiased individual 
interdependence could be of use to tackle problems such as how the use of an accessory in a group will be affected by the 
use of another when it is considered to be fashionable or unfashionable to wear them together. 

The study of the unbiased case is also useful to get some insight on what will happen in a more general situation 
of biased homogeneous groups. The introduction of nonzero constant opinion fields, will in general determine the sign 
of the average magnetisation and break the equilibria degeneracy, and the unpolarised state will no longer be stable. 
What were unpolarised regions for the unbiased case, will become characterised by an equilibrium demand or acceptance 
univocally determined by the constant IWA (or IWB and price). In the already polarised regions before the introduction 
of fields, however, multiple possible equilibria and first order phase transitions can be expected at low enough IWAs and 
temperatures [57] , 

Figures [1] to [3] show numerically calculated solutions which are stable (minima of the free energy). Green dots are used 
for the paramagnetic or unpolarised phase (only one stable equilibria) , blue asps (x) for ferromagnetic or polarised phases 
with both average choices of the same sign and crosses (+) for ferromagnetic or polarised phases with average choices of 
opposed signs (two physically equivalent stable equilibria in each case). Red triangles show mixed phase solutions (also 
two equilibria). When both types of ferromagnetic or polarised phases are superposed, there are four different equilibria, 
two stable (those aligned according to the sign of k) and two metastablc. 

The mixed segment is shown with a dashed black line (figures [5] and 13]). This is the only region where mixed phases 
are solutions and is the same for both models as these can only exist for k = (uncoupled case) and when one of the 
Ising models is in the paramagnetic phase while the other is still in the ferromagnetic phase, i.e., Jt < fi^ 1 < J s when 
Jt < J s (or J s < < J t when J s < J t ). Whenever there is a nonzero value of k, it is no longer possible for the opinion 
to remain unpolarised if the other group/opinion to which it is coupled is not. 

For the group interdependence or nonlocal model, the degenerate curve k 2 = J s Jt is drawn as a solid thick black line 
(figures [T] (a), [5] (a) and [3] (a) and (c)). As shown in the figures, there are no stable solutions of any type in the strong 
coupling regime. Frustration prevents the system from arriving at equilibrium. This was already being indicated by 
the fact that neither the completely polarised (s — ±I,t = ±1) solutions at zero temperature or the unpolarised one at 
very large temperatures are stable (which can be checked analytically). When considering that the interaction between 
both groups is of pure social type, assuming stronger social pressure from within the group than from the outside seems 
pretty natural. It can be considered a nice feature of the model that it breaks down when the problem is, sociologically 
speaking, ill defined (groups have not been chosen appropriately). The coupling between both groups, however, needs 
not to be of social nature. In this case, this would be pointing at a real impossibility of arriving at a stable equilibria if 
the interdependence between both groups is the main driver of the choice making process. 
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Figure 1: J t — (3 1 section for the (a) nonlocal and (b) local models for fixed J s = 1 and k = 0.3. 



The critical curve, where a second order phase transition takes place, is drawn in thin solid black (figures [TJ [5] and [3] 
(a) and (b)). This is given by the condition > max{J s , J t } together with (|2.12|) for the nonlocal model or (|3.9|) for the 
local model. It separates regions of polarised opinion (ferromagnetic phases) from others where the paramagnetic phase 
is the only stable equilibria (half of the population deciding each way for both choices for sufficiently large statistical 
fluctuations). When nonzero, constant IWAs are introduced, this will still separate regions where there will be a single 
possible value for the average choice from those where metastable equilibria can be expected to appear for low fields. 
Under such circumstances, drastic and irreversible shifts in the average opinion can be expected |57j . 

Figure [T] shows a J t — /3 _1 section for the nonlocal (a) and local (b) models for J s = 1 and k = 0.3. They arc 
remarkably similar qualitatively and quantitatively. If both groups/choices were not coupled, the critical curve would be 
the line of slope one across the origin and in all of the ferromagnetic region there would be two possible values (of same 
absolute value) for each average magnetisation (four possible states when combined), all with the same probability. The 
introduction of interdependence between both groups/choices promotes polarisation, making larger statistical fluctuations 
(temperature) needed as compared to the uncoupled case for imitation to promote a certain collective trend. For high 
enough social permeabilities /?, there will no longer be a threshold value of the social influence bellow which the system 
is in a paramagnetic phase. In the individual dependence case, even when we consider no imitation whatsoever for one of 
the choices (Jt = 0), social pressure on the other is enough for polarisation to appear. This is not the case of the group 
dependence model, where zero intra-couplings arc forbidden, as they always lay in the strong coupling regime. 

When considering ferromagnetic regions, the introduction of the coupling will, in general, fix the relative signs between 
both choices, although not the direction. For example, if we consider the adoption of a certain technical term by two 
professional groups or academic disciplines, if these are close and there is affinity between them, the new term will either 
gain acceptance or not in both groups. Similarly, if both groups perceive each other as confronted in any way, a high 
acceptance in one of them will always come with a low acceptance in the other. Likewise, if we analyse teenage outfit, if 
two accessories are perceived as fitting well, they will either both be used by over half of the population or by less than 
half of it, with no possibility of one of the trends being successful without the other. 
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For low enough temperatures and/or strong enough social pressure (coupling J t ), there is a region of metastability 
and the possibility of hysteresis. This behaviour is related to a first order phase transition at k = 0. If the inter-coupling 
is slowly reversing its sign, the system may not immediately change in the alignment of average choices, i.e., the system 
may get trapped in a local minima of the free energy (for a while, statistical fluctuations will make the system decay to the 
real ground state eventually). This is never the case for uncoupled groups/choices in the unbiased case and can be given 
interesting interpretations when considering two groups/choices with mutual influence slowly changing from negative to 
positive (or positive to negative). In the example of vocabulary adoption, this will imply the possibility of the use of the 
word remaining low in one of the professional sectors even as the mutual perception becomes positive, before there is an 
abrupt change and becomes widely accepted or not in both groups. In this region, the introduction of weak nonzero fields 
will probably enrich the metastability landscape |57j . 

Note that for the nonlocal model the metastability region extends up to /3 _1 =0 but it does not in the local case (which 
can be easily checked analytically). This difference can also be explained naturally in the light of social applications. Zero 
temperature is equivalent to a completely deterministic picture. Then, in the local case, each rational, utility maximiser 
knows exactly whether it is better off aligning or not both of its choices (although not in which direction) . In the nonlocal 
case, however, although there is still a globally preferred condition of alignment or not between both groups, the choice 
is not up to every agent as the coupling depends on what the other group is doing. 

Figure [2] shows k — sections for both models when J s = 1 and Jt = 0.6. The uncoupled case is described by the 
k = axis. Besides the promotion of polarisation and the differences already noted when examining the J t — /3 _1 section, 
for the local model, the critical curve is such that Wm^^^ fi^ 1 = J s + J t = 1.6. No similar asymptotic behaviour is 
shown by the nonlocal critical curve (for which high values of \k\ are of little interest as there are no stable solutions 
in the strong coupling regime). When the interrelation between both decisions is much more important than any of the 
imitation strengths, the critical value of social permeability at which a different qualitative picture can emerge due to 
social interaction, is basically determined by the intra-couplings. Further increasing \k\ will have little impact in the 
critical value of the temperature and, once permeability is fixed and while in a high \k\ range, small variations in the 
interdependence of both choices will never bring about a qualitative change in the average outcome. 

For both models, if /3 1 > J s + J t , the only stable state, if any, is the unpolarised one. In the nonlocal or group 
interdependence model, this is due to the intersection of the degenerate curve and the critical curve, precisely at J s + J t . In 
the local or individual interdependence model, it has to do with the asymptotic behaviour of the critical curve described 
above. This also makes sense in our context: at low enough permeabilities, in the local model, the interdependence 
between both choices can be made as large as desired without making social influence relevant qualitatively. This means 
that once we consider constant, nonzero IWAs, if the statistical fluctuations considered are large enough (such that 
(3- 1 > J s + J t ), both average choices will be aligned with the fields, no matter what the value of k is, and no multiple 
equilibria or abrupt changes will be possible. 

As for Jt — k sections, the situation is qualitatively different depending on the temperature chosen. Figures [3] (a) and 
(b) show the sections for the nonlocal and local models when = 1.5 > J s = 1. In this case we have mixed segments 
for k = and Jt > /3 _1 = 1.5, paramagnetic phase (unpolarised collective opinions) bellow the critical curve and no 
metastability region. As before, the uncoupled case is described by the k = axis. Again, only the local case presents 
asymptotic behaviour of its critical curve limifei^oo Jt — > /3 _1 — J s = 0.4. As there, for J t < — J s , only the unpolarised 
state will be stable (if any) for both models. If the dependence between both choices k is strong enough, the critical value 
of Jt will be practically independent of its value. 

Figures [3] (c) and (d) show the sections for the nonlocal and local models when = 0.4 < J s = 1. Mixed segments 
lie, for k = (uncoupled case), on < Jt < /3 _1 = 0.4. There is no paramagnetic region (both opinions will always 
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Figure 2: k — fi~ x section for the (a) nonlocal and (b) local models for fixed J s = 1 and J t = 0.6. 

be polarised). There is a region of metastability present for low enough |fc| and/or high enough J t , with the same 
socioeconomic implications described above. 

5 Conclusions 

Two different coupled Ising models have been analysed in the mean field regime in the context of binary choices and 
socioeconomic interactions. One of them is naturally suited to study how two groups affect each other when both arc 
making the same choice (group interdependence or nonlocal model). The other to the case of a group where individuals 
make two choices that are interrelated (individual interdependence or local model). Phase diagrams have been discussed 
for the zero choice fields case (unbiased populations). These can be useful when analysing the success of certain trends 
or traditions and is enough to delimit regions of the phase diagram where social and choice interaction will be important 
to the qualitative outcome when the fields are turned on. 

When compared to the uncoupled unbiased case, the introduction of interdependence, to an exogenous average choice 
or to another decision the individual is making, shifts the system towards polarisation, that is, larger statistical fluctuations 
(or smaller social permeability) is required for the average choices to remain unaffected by the social and choice interaction. 
Either both groups/choices will be simultaneously polarised or unpolarised. In the former case, there will in general be, 
as in the uncoupled case, two possible values of opposite sign for each average choice. The difference is that now the 
relative sign between them is in general fixed by the introduction of the coupling. This bonds the success of any fashion 
or trend to the success (positive interdependence) or failure (negative interdependence) in the other group/fashion to 
which it is coupled. There is a region at low inter-coupling and fluctuations (physical temperature) where it is possible 
for both average choices to align differently from the coupling between them. This allows for sudden, irreversible shifts 
from low to high (or high to low) demand or acceptance (which was never the case in the uncoupled unbiased case), and 
for strong history dependence. 

When comparing the effects of coupling the choice to the same decision making process in an external group to that of 
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another choice each individual is making, there are also some interesting differences. The coupling of groups will prevent 
them from arriving at a state of stable equilibrium whenever it is the main driver of the decision making process. In 
the completely deterministic picture, only group interdependence will allow for the possibility of alignment between the 
average choices contradicting the sign of the coupling. Although, contrary to the case of group interdependence, there 
are stable states for high individual choice interdependence, at large enough values of the inter-coupling, the polarisation 
threshold will be independent of it and determined by the sum of social imitation on each choice. 

The models presented are a first look into the problem and share the limitations exposed for the uncoupled case 
in the introduction. Probably the greatest weakness of this approach is that of considering completely homogeneous 
populations. The first natural way to encode a group's particularities and preferences seems to be the use of random 
choice fields, which is part of work in progress [57]. From a wider perspective, the most important line of work might be 
testing the suitability of these and other related models to explain -ideally quantitatively but at least qualitatively- data 
for real social processes. This is no simple task and there is currently also work being done on these lines |58j . 
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